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On the nature of an ejection event in the jet of 3C 111 
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Abstract. We present a possible scenario for the ejection of a superluminal component in the jet of the Broad 
Line Radio Galaxy 3C 111 in early 1996. VLBI observations at 15 GHz discovered the presence of two jet features 
on scales smaller than one parsec. The first component evolves downstream, whereas the second one fades out 
after 1 parsec. We propose the injection of a perturbation of dense material followed by a decrease in the injection 
rate of material in the jet as a plausible explanation. This scenario is supported by ID relativistic hydrodynamic 
and emission simulations. The perturbation is modeled as an increase in the jet density, without modifying the 
original Lorentz factor in the initial conditions. We show that an increase of the Lorentz factor in the material 
of the perturbation fails to reproduce the observed evolution of this flare. We are able to estimate the lifetime of 
the ejection event in 3C 111 to be 36 ± 7 days. 
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1. Introduction 

Flaring events at radio frequencies are known to take place 
in Active Galactic Nuclei (AGN), usually followed by the 
obser vation of new radio fea tures in the parsec-scale jets 
(e.g., Savolainen et al. 2002[ ). It has been shown that the 
ejection of those features, or components, is related to 
dips in the X-ra y emission from the a ctive nucleus in the 
case of 3C 120 (IMarscher et al.ll2002l ). and perhaps also 



ciated with higher radio brightness, to the soft/high state, 
associated with lower radio emission - a decrease in the 



inject ion of emitting particles in the jet (jFender fc Belloni 



2004 and references therein). 



in 3C 111 (jMarscheij l2006a[ ) . The dips in X-rays precede 
the observations of new radio-components. The decrease 
in X-ray emission may be caused by the loss of the in- 
ner regions of the disc. In this scenario, a fraction of the 
accreted material is injected in the jet and a new com- 
ponent is later observed in VLBI images, after the ma- 
terial becomes detectable at the observing frequencies, as 
it evolves downstream. The components are interpreted 
as the shocks produced by the ejection of denser and/or 
faster plasma in the flar ing event from the accretion disc 
( Marscher fc GeajllQSSTl . The conditions for triggering the 
ejection of the material in those radio features are still un- 
known. In the case of microquasars, it has been proposed 
that the stronger components are ejected right before the 
passage of the source from the X-ray hard/low state, asso- 
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If we interpret jet components as shocks propagat- 
ing in a supersonic flow, their origin must be related 
to an increase of pressure and/or velocity in an in- 
jected perturbation with resp ect to the steady initial flow. 
Hydrodynamical simulations ( Alov et al]|2003l AOS here- 
after) have shown that such jet perturbations produce a 
forward and a reverse structure, which would be expected 
to be observed as a fast front and a slower back compo- 
nent. 

In the jet of the Broad Line Radio Galaxy (BLRG) 
3C 111 {z = 0.049, Imas ~ Ipc), a very strong flar- 
ing event in early 1996 gave rise to the ejection of two 
jet features observed at 15 GHz with the Very Long 
Baseline Array (labeled as compone nts E and F - see 
Fig.fDand lKadler et al.ll2008L hereafterlKOS') . Both compo- 
nent trajectories can be back-extrapolated to similar ejec- 
tion epochs within 3 months (around 1996.10). However, 
they show different speeds and the time evolution of their 
brightness is different (see Fig. [T]): the inner component 
F is initially brighter (1996.82 and 1997.19) and fades 
out very rapidly (1997.66 and 1998.18), while the lead- 
ing component E shows a slower decrease in flux den- 
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Fig. 1. Core distance and flux density evolution with time 
of co mpon ents E and F in 3C 111, based on the results 
fromlKOSl. 



sity. After 1999 (see KOSi ). F disappeared and E evolves, 
accele rating and genera ting trailing components in its 
wake ( Agudo et al.l2001 ). The differences betweeen E and 
F cannot be attributed to different Doppler factors of 
the components, a s the se are very similar {De ~ 3.2, 
Dp '-^3.1, foUowingiKOS)- The possibility that component 
F represents a second injection after component E is highly 
improbable on the basis of its velocity and its brightness 
evolution, which cannot be linked to the propagation in 
the wake of the latter. In this letter, we investigate, in 
a qualitative way, the possibility that these components 
are the front and rear region of a single perturbation. A 
quantitative comparison between the simulations and the 
observations would require detailed knowledge about the 
nature of the flow and is out of the scope of this work. 
The numbers and error estimates used in this letter that 
are related to the observations are taken or derived from 



K08 



2. Hydrodynamics and emission 

In numerical simulations, the enhanced injection of mate- 
rial in jets has been modelled as an abrupt square per- 
turbation of the flow d ensity and/or Lorentz factor at 



the in jection point (e.g., iGomez et al.lll997t lAgudo et al 



2001 ^ in a steady jet that recovers the initial state after 
the perturbation is completely ejected. A03 show that the 
perturbation generated by the injection of denser plasma 
propagates downstream, spreading in the axial direction 
along the jet and finally splitting into two distinct regions. 
The forward region moves against the underlying flow and 
the reverse region propagates backwards in the jet fluid 
reference frame, thus acting against the fluid that is in- 
jected after the perturbation. The injection of dense fluid 
behind the perturbation feeds the reverse shock, making 
it potentially observable (with similar brightness to that 
of the forward shock, A03) for relatively long times. 

We have performed one-dimensional numerical rel- 
ativistic hydrodynamic (RHD) simulations in which a 
square perturbation in density is injected into a steady 
jet, without modifying the initial Lorentz factor, and re- 
laxing the condition that the initial jet flow is reestab- 
lished immediately after the perturbation. We have sub- 
stituted this by a rarefied flow, representing a reduction 
of the injection rate, in order to avoid the formation of a 



strong reverse shock. In this picture, the original jet injec- 
tion rates should be recovered after some time. However, 
in this work we only focus on the evolution of the strong 
ejection and the period before the reestablishment of the 
jet flow. Multidimensional simulations are out of the scope 
of this work due to the computational effort required and 
to the one-dimensional character of this problem. The sim- 
ulations have been performed using a numerical code that 
solves the equations of relativistic hydrody namics written 
in the conservation form, as described in I Perucho et al 



(|2005l ) and iMarti et al.l (|1997^ . The details of the simu- 
lation are given in the caption of Fig. ^ The top panels 
in Fig. [2] show different snapshots of the evolution of the 
square perturbation injected in a steady flow, in pressure, 
Lorentz factor and specific internal energy. Using the RHD 
simulations as input, we have computed the corresponding 
ID optically thin radio synchrotron emission as seen by a n 
observer with a line of sight at 19° to the jet axis ( KOSi ). 
For these computations, w e used the numerica l code and 
the procedure described in lGomez et al. ( 1997 ) and refer- 
ences therein. This code takes into account all the relevant 
relativistic effects, including the light travel time delays. 

In the simulation (see Fig. [2]), the front region includes 
the leading part of th e perturbation and is identified with 
component E in whereas we define the fading region 
as the rear part of the perturbation and identify it with 
component F (see Fig. [5]). The material in the front region, 
consisting of shocked material from the steady jet and rar- 
efied material from the perturbation separated by a con- 
tact discontinuity, shows smaller values for the pressure, 
and some acceleration due the propagation in the lower 
pressure steady jet fluid. The material in the fading region 
crosses the receding rarefaction that separates it from the 
front region (top panels in Fig. [2]) and it is also "eroded" 
by the back rarefaction. Consequently, the front structure 
evolves, increasing its size as the front shock incorporates 
material from the steady jet and the material from the fad- 
ing region crosses the receding rarefaction. Thus, the front 
region consists of the forward shock structure of the per- 
turbation (E in Fig. [2|) , and the fading region is formed by 
the remains of the perturbation that have not crossed the 
receding rarefaction (F). The synchrotron emissivity (bot- 
tom panel in Fig. [5]) is governed by the jet pressure and 
hence the emission evolution is very similar to the pressure 
evolution of the RHD simulations. In the emission results, 
the front region (component E) propagates without much 
flux density evolution after injection. However, the fading 
structure (component F), which initially shows a notably 
larger flux density than component E, rapidly decreases 
in emission as the receding and back rarefactions erode it. 
The reverse shock (see A03) is neither relevant nor obser- 
vationally significant in our simulations, as it propagates 
in a very rarefied medium. For this reason it is not shown 
in Fig. m 

Notice that the Lorentz factor values in Fig. [2] are 
those corresponding to the fluid. In contrast, VLBI ob- 
servations provide us with pattern velocities. In the sim- 
ulation, the velocity of the front shock is measured to be 
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Fig. 2. Snapshots of the evolution (left to right) of a square perturbation injected in a steady jet, followed by a 
strong rarefaction. The dotted-light-blue lines stand for the Lorentz factor, the solid-dark- blue line for pressure and 
the dashed-red lines for specific internal energy. The simulation is run with 24000 cells; the velocity of the initial flow 
is Vj = 0.9 c; the perturbation is injected during a time interval 0.2 Rj/c, with a density twice that of the jet, the same 
specific internal energy and velocity Vp — 0.9 c; the rarefied medium is injected after the perturbation with the same 
velocity as the initial flow, and a pressure ten times smaller than that of the initial flow. Please note the change of scale 
in the abcissae. The bottom panels show the simulated total intensity emission along the jet axis at fo ur rep resentative 
epochs. The identification of the features in the simulation with the observed components E and F in IkosI is indicated 
in each panel. A jet width of 500 cells and axial symmetry is used to compute the emission. 



Vs ^ 0.96 c {v^^ ~ 3.5c), whereas that of the fading region 
is Vr ~ 0.87 c (wp*^ ^ 1. 7c), both similar to those found 
in the observations ( KOSi ). The velocity of the material in 
the fading region is faster than that of the receding rar- 
efaction (cf. Figl2]), as expected from the explanation in 
the previous paragraph. We also point out that the dilute 
material shown in Fig. [2] presents a modified velocity due 
to passage through the reverse shock. 

A second simulation was performed for a faster per- 
turbation, with Lorentz factor F = 3.6, while keeping the 
rest of the parameters as in the previous simulation. The 
results (see Fig. [3]) show that the front region of the per- 
turbation is overpressured with respect to the rear region 
and, thus, the former is brighter than the latter, as shown 
by the emission simulations (bottom panels in Fig.[3|). This 
is in clear contrad iction to the observations of the jet in 
3C 111 (Fig.[T]and|K0i). The difference is due to the pres- 
ence of a stronger front shock. It is also important that 
the wave separating both regions is now a reverse shock, 
instead of the receding rarefaction shown in Fig. [2] This is 
a general result for fast perturbations, including the case 
of a fast perturbation in pressure equilibrium with the 
steady jet. 

In any of the scenarios given above, the perturbed re- 
gions have enhanced emission with respect to the under- 
lying jet. However, only in the case of an overpressured 
perturbation with the same Lorentz factor as the under- 
lying flow, and including the presence of a rarefaction be- 
hind the perturbation, does the second component fade 
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Fig. 3. Same as Fig.[2]but showing the evolution of a fast 
square perturbation injected in a steady jet followed by a 
strong rarefaction (top panels) and representative emis- 
sion plots (bottom panels). The perturbation is injected 
with a density twice that of the jet and Lorentz factor 
F = 3.6. This image compares with the first two snap- 
shots in Fig. [2] and shows that, in this case, component F 
would not be brighter than component E. The later stages 
of evolution in this simulation are very similar to those in 
Fig. [2] and are not shown here. 
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out rapidly and then the first one dominates the emission, 
as observed for components E and F in 3C111. 

3. Discussion and conclusions 

We are able to describe the evolution of the leading com- 
ponents of a major ejection in the BLRG 3C 111 as a 
perturbation of dense (overpressured) material followed 
by a dilute medium. The qualitative explanation given 
in the previous section supports the picture of the ejec- 
tion of perturbations of dense material giving rise to radio 
components. Our model would remain valid if the ampli- 
tude or time-duration of the perturbation is changed. In 
addition, we postulate that these ejections could be fol- 
lowed by a decrease in the injection rate of bulk flow par- 
ticles. This avoids the formation of a reverse shock, which 
would lead to a qualitatively different observational result. 
Observational support for the inclusion of this tenuous 
material in the simulations can be found in the promi- 
nent emissio n gap following behind the E/F complex in 
3C 111 fsee iKOSf ). New ejection of emitting material is 
detected on the time scale of more than 2 years, corre- 
sponding to a gap width of up to 2 mas in 1999 (cf. IK08 ). 
The observations of microquasar jets show that, in gen- 
eral, major ejections are followed by a decrease in radio 
brightness ( Fender fc Bellonl 2004 ). If wc interpret this as 
a decrease in the injection rate in the jet, while the sys- 
tem relaxes back to the initial steady state or generates the 
conditions for the injection of a new perturbation, it could 
be a process similar to that explained here. However, the 
strong decrease in radio emission observed in microquasar 
jets is not observed in AGN jets. This setup, in which we 
put in relation the processes taking place in the jet and 
the accretion disk, relies also on th e results of multiwave - 
len gth observational c ampaigns by lMarscher et aL ( 2002[ ) 
and lMarscheii (j2006al ). The latter work showed a relation- 
ship between dips in the X-ray emission from the accretion 
disk and the ejection of radio components in 3C 120, and 
possibly in 3C 111. 

We can place an upper limit for the lifetime of the ejec- 
tion event in its passage through the radio core (opaque 
and compact emitting region at the origin of the radio 
jet) using only observational data. We consider that: a) 
the perturbation starte d to cross the radio core in the ob- 
served jet in 1996.10 (|K08[ ). b) the receding rarefaction 
"eroding" component F moves with the observed veloc- 
ity vp = 0.91c (Figd]), and c) the last epoch at which 
component F is observed (1998.18) is taken as the time 
at which the receding rarefaction has completely eroded 
this component - which is justified since the flux of com- 
ponent F is one order of magnitude smaller than that of 
component E in this epoch (K08) . With these assump- 
tions we can calculate the time that the last material of 
the perturbation needs to catch up the receding rarefac- 
tion, and from that, we can estimate the time lapse of the 
crossing through the core. We use the velocity of compo- 
nent E as an upper limit for the velocity of the fluid. The 
result tells us that the crossing of gas through the core 



had to end as soon as (1.98 ± 0.02) years before epoch 
1998.18, this is, (0.10 ± 0.02) years after 1996.10. This 
means At ^ (36 ± 7) days between the passage of the 
first material and that of the last portion of gas in the 
perturbation. This lapse of time should be smaller if the 
gas catching up the rarefaction is slower than the velocity 
of the head. This represents a first order estimate, which 
only depends on the ratio between the velocity of the ma- 
terial in the perturbation and that of the rarefaction wave, 
and provides an upper limit for the duration of such an 
event. We could consider this lifetime of the ejection as an 
upper limit for its triggering event in the black- hole/inner- 
accretion-disc system, as we are not taking into account 
the coUimation and acceleration processes, which should 
be kinematically important in the most compact scales. 
In the case of microquasar jets, the radio flares that fol- 
low the dips in X-rays -and are related to the ejection 
of components- have been e stimated to last between sec - 
onds and tenths of minutes ( Mirabel fc Rodriguez Il999f ). 
i.e., lO"' — 10^ times smaller time-scales than for 3C 111, 
whe reas the ratio of the black-hole masses is of the o rder of 
10^ (iMirabel fc Rodriguejll999l : lGrandi et"aDl2006l) . If we 
take into account that our result represents an upper limit 
for the timescale of the ejection event, this suggests that 
the time-scale factor between quasars and microquasars 
may not come directly from the ratio of black hole masses 
(see also lMarscheiil2006bl ) . 

Strong radio flares associated with the ejection of 
radio-components should be carefully followed up, with 
a sufficiently dense time sampling, in 3C 111 and other 
sources for which we can achieve similar or better lin- 
ear resolution with the VLBI technique. This would help 
in performing analyses like that presented here and to 
test our conclusions. The stretching of the size of any 
relativistic structure propagating through a jet must 
largely favour the detection of such double structures in 
the jets of nearby AGN. Databases provided by moni- 
toring programme s such as the MOJAVE/2 cm VLB A 
([Kellermann et al.l l2004t iLister fc HomanI 120051 ) survey 
are very valuable in this context. Future work in- 
cludes monitoring in radio/X-ray campaigns of different 
AGN sources and further numerical calculations includ- 
ing multi-dimensional RHD, RMHD and emission simula- 
tions. 
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